Characteristics of the asymmetric simple exclusion process in the presence of 

quenched spatial disorder 
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We investigate the effect of quenched spatial disordered hopping rates on the characteristics of the 
asymmetric simple exclusion process (ASEP) with open boundaries both numerically and by exten- 
sive simulations. Disorder averages of the bulk density and current are obtained in terms of various 
input and output rates. We study the binary and uniform distributions of disorder. It is verified 
that the effect of spatial inhomogeneity is generically to enlarge the size of the maximal current 
phase. This is in accordance with the mean field results obtained by Harris and Stinchcombe [22l ]. 
Furthermore, we obtain the dependence of the current and the bulk density on the characteristics of 
the disorder distribution function. It is shown that the impact of disorder crucially depends on the 
particle input and out rates. In some situations, disorder can constructively enhance the current. 



I. INTRODUCTION 



Transport processes in disordered media constitute an 
important class of problems especially in the light of 
their relevance to the modelling of a vast variety of phe- 
nomena in physics and many interdisciplinary areas. A 
partial list of applications includes transport phenomena 
in porous media, diffusion in biological tissues and con- 
duction through composite solids P, 0, Q ■ It is a well- 
established fact the disorder can strongly affect the trans- 
port characteristics of equilibrium as well as out of equi- 
librium systems. Among various non equilibrium sys- 
tems, low dimensional driven lattice gases have played 
an important role in describing the transport in many 
physical, chemical and biological processes [J, d, 0, 0, Hf . 
In particular, one dimensional driven diffusive systems 
in the absence of disorder have been extensively studied 
during the past two decades and at present there exists 
a rich literature of results both analytic and numeric Q . 
Phase structures of these systems are well known. It is 
well understood that non-equilibrium systems can exhibit 
phase transitions in low dimensions. A model which has 
played a paradigmatic role in out of equilibrium statis- 
tical physics is the asymmetric simple exclusion process 
(ASEP) @. The model is amenable to exact analyti- 
cal solution [l(| EH, EH- Therefore it is a natural and 
important question to investigate the effect of quenched 
disorder on the phase structure of ASEP. Recently some 
efforts and new strides have been made in the challenge 
between disorder, interaction and drive. The exploration 
of the disordered ASEP began with a single defective site 
in a periodic chain by Janowsky and Lebowitz (T3 , [lij . 
They showed that even one defective site can remarkably 
lead to global effects on the system current and density 
profile. Evans solved the ASEP with moving impurities 
where particle hopping rates were chosen randomly from 
a distribution function 15]. It was shown that special 
distribution functions can give rise to a new phase tran- 



sition analogous to Bose condensation. Subsequently, 
Tripathy and Barma [H, E3 considered the ASEP on 
a ring with many defective sites. Their investigation re- 
vealed the existence of phase segregation in a wide range 
of global densities in the chain. In conjunction with the 
results of the ASEP on a ring, an investigation of the 
disordered ASEP in an open chain was introduced by 
Kolomeisky [l8|. He showed that in some ranges of in- 
put and output rates, a single defect in the bulk could 
affect the system properties on a global scale. Recently a 
new wave of attention has been created on the disordered 
ASEP[H H3, EH, EH, SI H 11, HE M, H|. In partic- 
ular, Chou and Lakatos have studied the effect of a few 
defective sites in the open ASEP (23[. Their investiga- 
tions have revealed that generically the disorder's impact 
is highest when the number of defects is very small. In- 
creasing the number of defects above a certain value has 
no further effect on the system current. The question of 
the effect of a single defect in the ASEP coupled with 
a 3D bulk reservoir with adsorption/desorption kinetics 
was recently addressed by Frey (29|. Besides, some time- 
dependent aspects of the disordered ASEP has been dis- 
cussed by Barma [301 ] . Our goal in this paper is to deal 
in some more depth with the problem of the disordered 
ASEP. Especially we will focus on the role of the binary 
distribution function where it provides the possibility of 
simultaneous study of both the strength and the den- 
sity of disorder throughout the chain. Via extensive MC 
simulations, we show for the binary and uniform distribu- 
tion functions the generic impact of disorder is to reduce 
the size of the low and high density phases. More in- 
terestingly, we show in some circumstances, disorder can 
constructively act in a manner to increase the system 
current. 
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FIG. 1: Asymmetric exclusion process with spatial disorder. 

II. DESCRIPTION OF THE PROBLEM AND 
NUMERICAL SOLUTION 

To keep the paper self-contained, let us first define the 
disordered ASEP. fmagine a one dimensional stochastic 
process defined on a discrete ID lattice of length L. Each 
site can hold at most one particle. We assign an integer 
valued number Si to each site i (see figure (I ) ). If site i is 
occupied, Sj=l. If it is empty then is zero. The system 
configuration at each time t is characterized by specifying 
the occupation numbers Sj i = 1, ■ ■ ■ , L. During an 
infinitesimal time dt each particle can stochastically hop 
to its rightmost neighbouring site provided the target 
site is empty. If the target site is already occupied by 
another particle, the attempted movement is rejected. 
The hopping takes place with a site dependent rate pi 
which is drawn from a given distribution function f(p). 

There is no spatial correlation between the set pi i — 
1, ■ ■ ■ , L and correspondingly the p^s can be regarded 
as independent stochastic variables which are identically 
distributed according to the site independent distribu- 
tion function f(p). Denoting the averaged local density 
at site i by n,, one can simply write the following rate 
equations for a particular realisation of hopping rates 
Pi,- ■■ ,Pl (i = 2, •• • ,L- 1): 

=p i -i(n i - 1 (l - rii)) -pi{m(l - n i+ i)), (I) 



-(m) = a(l - (m» - Pl (ni(l - n 2 )), (2) 

^(n L ) = p L -i{n L -i(l - n L )) - (3{n L ). (3) 

No exact analytical solution exists for the above set 
of nonlinear differential equations. Restricting ourselves 
to the stationary state properties of the system, we set 
the left hand sides equal to zero. We further simplify 
the equations by taking the assumption of a mean field 
equation where the two-point functions are replaced by 
the product of two one-point functions. This assumption 
reduces the steady state equations into a set of nonlinear 
algebraic equation with L unknowns m, 712, ■ • • , n L . 

A. Numerical approach to mean field equations 

Even by employing the assumption of a mean field, 
we are not able to solve the nonlinear algebraic equa- 
tions. Therefore, we should resort to numerical methods. 



We now outline a numerical approach for solving the set 
of nonlinear equations. The approach is based on the 
shooting method for solving boundary value problems. 
To this end, we choose a trial n\ denoted by n* r and 
successively evaluate n 2 , • • • ,Ul through forward itera- 
tion. The system current then turns out to be a(l — n* r ). 
Since the current should be equal for all sites, if the 
guessed value of was correct then the current eval- 
uated from the last site i.e., /3n^ would have the same 
amount a(l — n\ r ) evaluated from the first site. To 
match these currents, we gradually increase n l { from zero 
and evaluate both currents from the first and last sites. 
Whenever these two values become equal, then we have 
a solution. Note that in an acceptable solution, all the 
densities (m), • • ■ , (ul) should lie between and 1. We 
are interested in knowing the overall effect of disorder 
on the transport characteristics of the ASEP. For given 
values of a and /?, we evaluate the current and density 
for many samples of disordered chains and average over 
these samples. We denote the sample averaged current 
and bulk density by (J) and (p) respectively. For ob- 
taining a better insight, we have also executed extensive 
Monte Carlo simulations. The disorder distributions we 
consider consist of uniform and binary. More explicitly, 
the normalized uniform distribution in the interval [a, b] 
has the functional form fi(p) = 53^ with mean and 

variance ^ b ~$ ■ The binary distribution has the form 
h(p) = fd(p—pi) + (l — f)5(p—p2) where the binary rates 
Pi and pi and their probabilities / and 1 — / are given. 
The mean value and the variance are fp\ + (1 — f)p2 and 
fPi + (1 - f)p\ - [fPi + (! - f)P2? respectively. In the 
subsequent sections, we show the result of simulation as 
well as numerical solution of the mean field equations. 

III. BINARY DISTRIBUTION FUNCTION 

Let us first discuss the binary distribution of the 
quenched disorder. Although in this type of distribu- 
tion the defect strength is allowed to take only an inte- 
ger number of values (here 2), but even in this simple 
case one encounters some nontrivial aspects which are 
worth investigating. Figure (2) depicts the dependence 
of average current versus (3 for some fixed a. The pa- 
rameters of the binary distribution function is as follows: 
/ = 0.5, p\ = 0.8 and P2 = 1 .2. The mean value of the 
quenched hopping rate is fixed at I. The bulk density 
and the current have been averaged over 1000 disordered 
samples and the system size is 200. 

One observes similar behaviour to the normal ASEP. 
Currents rise up to a critical (3 and then get saturated. 
The overall effect of disorder is to reduce the value of 
the currents in each phase. In the normal ASEP, the 
dependence of current on (3 in the high density (HD) 
phase is (3(1 — (3). Saturation of current means that we 
are in the maximal current (MC) phase. However, the 
current saturates at 0.15 which is less than the value 
of the maximal current 0.25 in the normal ASEP. The 
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FIG. 4: Disorder-averaged densities vs for various input 
FIG. 2: (J) versus f3 for various input rates a. ra tes a. System size is 200 and averaging has been executed 

over 1000 disordered samples. 
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FIG. 3: Disorder- averaged current vs a for various output 
rates (3. System size is 200. 



reason is due to the presence of defects which slow down 
the current. Figure (3) exhibits the dependence of (J) 
versus the input rate a. 

The behaviour seen in the above graph is analogous to 
the normal ASEP with the difference that the disorder 
has yielded to an overall diminishing of (J). In the nor- 
mal ASEP, the dependence of (J) on a in the low density 
phase is a(l — a). We note that upon entering the satu- 
ration regime i.e., the MC phase, the current value which 
is 0.15, is less than that of the normal ASEP 0.25. The 
dependence of bulk densities on a and f3 are exhibited in 
figures (4,5). Similar to current diagrams, the overall be- 
haviour is analogous to the normal ASEP. Here the effect 
of disorder is to enhance the densities. 

One observes the persistence of the first-order high to 
low density transition. The saturation density is slightly 
above the normal ASEP value a which is due to defects. 
For higher a corresponding to the MC phase, the limiting 
density is the same as the normal ASEP i.e., (p) = 0.5. 
This indicates that the presence of defects does not alter 
the density in the maximal current phase but reduces the 
current as discussed earlier. The dependence of (p) on a 
is shown in figure (5). 

Similar to fig. (4), the density value in the saturation 
regime is roughly 0.5 which is the same as the normal 
ASEP density value in the MC phase. However, (p) dif- 



FIG. 5: Disorder-averaged densities vs a for various output 
rates j3. System size is 200 and averaging has been executed 
over 1000 disordered samples. 



fers from 1 — j3 in the LD phase. We note that both high 
to low and low to high phase transitions which are first 
order in the normal ASEP are replaced with a smoother 
behaviour in the presence of disorder. We have exten- 
sively performed Monte Carlo simulations for all ranges 
of a and (3. The simulation results confirm the existence 
of three phases of LD, HD and maximal current. Further- 
more, our simulations show the growth of the maximal- 
current region and shortening of the sizes of the LD and 
HD phases. These findings are in agreement with the 
mean-field based conclusions of Harris and Stinchcombe 
■ If one changes the parameters of the binary distri- 
bution function, the overall picture remains qualitatively 
the same as in the above diagrams. Nevertheless, the 
quantitative values of both (J) and (p) in the phases 
depend on the parameters of the distribution functions. 
More concisely, currents and densities are functional of 
f(p). In the following figure, we exhibit the phase dia- 
gram of disordered ASEP for some binary distribution 
functions. 

We note that the size of the MC phase is an increasing 
function of the variance of the distribution function. The 
reason is that currents and densities are dominated by the 
number of defective sites. If the variance of the distribu- 
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FIG. 6: Phase diagram of the disordered ASEP for some bi- 
nary distributions. The distribution parameters are specified 
in the figure. All the distributions have / = 0.5 and (p) = I 
but the variances are different. 



tion is large, then the probability of finding sites with 
notably small hopping rates is considerable and there- 
fore (J) and (p) are highly affected. The enlargement of 
the MC phase has been reported for ASEP with a sin- 
gle defect in the bulk [l8j]. We recall from normal ASEP 
that the critical values of the input and output rates are 
a c = (3 C = | where p is the hopping rate. In principle, 
since a, (3 and p are rates, they can vary from zero to 
infinity. Therefore, it is possible to choose the time unit 
such that p scales to unity. In the disordered version of 
the ASEP, one does not have a single hopping rate so it 
would be better not to restrict ourselves to a particular 
time unit. For the sake of comparison we write the values 
of the density and current in the low density, high density 
and maximal current phases of the normal ASEP: 

a (3 1 , A . 

PhD = ~, PHD = 1 , PAW = r- (4) 

p p 2 

Having in the mind that J = pp(l — p) we thus obtain: 

ex. 



Jld = ol(1 ). 

P 



Jhd = /3(1 - -)• 
P 



Jmc - \- 



(•5) 



(6) 



(7) 



Since in the binary distribution there are three param- 
eters namely pi,P2 and /, we have studied two distin- 
guished cases. First, we restrict ourselves to the condi- 
tion (p) — fpi + (1 — ,f)p2 = 1- This leaves only two free 
parameters. In the second case, we impose the condition 
P2 = 1 while / and p\ are free to take arbitrary values 
in [0, 1]. In the latter case, the emphasis is on the role 
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FIG. 7: Figures 7-9: (J) vs / for various pv. a = 0.05, (3 = 
0.8 (top), a = 0.8, fi = 0.05 (middle) and a = 0.6, /3 = 0.8 
(bottom). System size is L — 300. 



of defective sites (pi < 1) among normal sites (p2 = 1) 
whereas the former case allows having fast hopping sites 
with rate p% > 1. Our simulations showed that there is 
no significant differences between the results of these two 
cases. Therefore in what follows, we only exhibit the re- 
sults for the case P2 — 1 i.e., slow defective sites among 
normal sites. Our first set of graphs (all obtained via 
Monte Carlo simulations) illustrates the dependence of 
(J) on / for various p\ in three sets of input and out- 
put rates corresponding to low input- high output, high 
input-low output and high input-high output. 

When a is small and (3 is high (fig. 7), the effect of 
increasing / is to reduce the current. Smaller values of 
Pi exhibit a sharper decrease. This is natural since the 
bulk density is low and therefore the system current is 
more sensitive to both the number and the strength of 
the defects. The dependence of (J) on / changes quali- 
tatively when one goes to the situation characterized by 
high a and low (3 ( fig. 8 ). Here we are confronted with 
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unexpected and novel features. For p\ less than 0.15, the 
current appears as a decreasing function of / while for 
Pi > 0.1 it increases up to a maximum and then starts di- 
minishing. Accordingly, the optimum value of / at which 
(J) is maximum is no longer / = but rather a nonzero 
fmax- This implies that the effect of disorder is to en- 
hance the current which is a desirable effect. The location 
of fmax shifts towards higher values when one increases 
P\. For large p±, (J) shows an increasing behaviour with 
a small slope. The slope tends to zero when pi — > 1. 
For relatively high values of a and (fig. 9), one still 
observes that the dependence of current versus / shows 
a decreasing character. In this case, the disorder has the 
expected behaviour i.e.; the higher the number of the 
impurities, the larger the decrease of (J). However, the 
interesting point is the abrupt change in the behaviour 
of the current reduction. For each p\, the current shows 
a rapid reduction up to a certain / and then decreases 
very smoothly in a nonlinear fashion. This marks the 
fact that impurities affect the system beyond a certain 
relative frequency. These results are in agreement with 
[28l |. To gain a deeper insight, it would be instructive to 
look at the behaviour of (p) vs /. In low input and high 
output rates (fig. 10), one observes that for low p\ the 
bulk density rises up to a maximum and then decreases 
even below the value of normal ASEP. One might naively 
think that increasing the density of defective sites leads 
to an enhancement of the bulk density due to the for- 
mation of high density regions behind them. However, 
the point is that if the density of defective sites reaches 
a certain value, the probability of finding defective sites 
in the vicinity of the first site of the chain increases too. 

This in turn gives rise to a blocking of the current of 
particles in the chain bulk. As a result, a large portion of 
the bulk remains almost in a LD regime which leads to a 
decreasing in the bulk density throughout the chain. This 
scenario remains valid for small p\. For larger pi, one 
observes the expected increase of the bulk density upon 
increasing the density of defective sites /. The reason is 
that once the defect strength is reduced below a thresh- 
old, the formation of high density regions behind these 
weak sites will be suppressed and therefore the particles 
can more easily flow throughout the bulk. As a result of 
this flow, enough particles can be found in the bulk. This 
increases the number of local high regions behind defects 
which in turn give rises to the enhancement of (p). In 
the high a-low regime (fig. 11) and for p\ < 0.2 the 
dependence of the bulk density on / is sharply decreasing 
up to a certain / . Afterwards, (p) becomes independent 
of / and a lengthy plateau region forms. At / = 0.9, (p) 
shows a rather linear increase to its asymptotic value 0.5 
in the MC phase. We note that when / = 1, all the sites 
are defective. For instance, in the case p\ = 0.1, the crit- 
ical input and output rates are a c = C — ^- = 0.05. In 
this case a = 0.8, = 0.05 lie in the MC phase and hence 
(p) approaches 0.5 in the limit / — > 1. The reason is that 
when the input is high and the output is low, impurities 
give rise to phase segregation behind them [13, • The 
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FIG. 8: Figures 10-12: Bulk density dependence on / for 
various values of pi: low density phase a — 0.05, (3 = 0.8 
(top), high density phase a = 0.8, f3 == 0.05 (middle) and 
maximal current phase a = 0.6, /3 = 0.8 (bottom). System 
size is L = 300. 



formation of macroscopic low density regions in front of 
them leads to a sharp reduction of (p). For p\ > 0.2 the 
decrease of (p) becomes much more smooth. The reason 
is that weaker defects are unable to produce low enough 
density regions. The other interesting point is that when 
the input rate is high, increasing the number of defects 
will prevent a high inflow of particles, which is due to 
the largeness of a and regulates the flow along the bulk. 
The overall effect is to reduce (p) from high a to much 
lower values. As we had already seen in the current dia- 
grams, in some values of p%, this diminishing in the bulk 
density is accompanied by the current increment as ex- 
hibited in fig. (8). Now let us discuss the regime where 
a and are both greater than 0.5 (fig. 12). For weak 
defect strength, the density is almost independent of /. 
For pi less than 0.5, (p) shows a smoothly decreasing 
dependence on / until it becomes independent of / and 
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FIG. 9: Figures 13-15: Current dependence on p\ for various 
values of /: low density phase a = 0.05, /3 = 0.8 (top), high 
density phase a — 0.8, /3 = 0.05 (middle) and maximal current 
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FIG. 10: Figures 16-18: (p) vs p\ for various /: LD phase 
a = 0.05, f3 = 0.8 (top), HD phase a = 0.8, f3 = 0.05 (middle) 
and MC phase a = 0.6, /3 = 0.8 (bottom). System size is 
L = 300. 



correspondingly a plateau region forms. The length of the 
plateau is relatively large and becomes larger for smaller 
values of p\. Increasing / beyond the plateau value, one 
again encounters an increasing behaviour of (p) until it 
reaches the normal ASEP value (p) = 0.5 in the limit 
/ = 1 in which all sites have become defective. In order 
to shed more light on our understanding, we now study 
the effect of varying the disorder strength p\ for fixed 
values of /. Analogous to the previous studies, we first 
consider the current which is shown in figures (13-15). 

When a is small and (3 large, the effect of increasing p\ 
is to increase the current to its normal value a(l— a). For 
each /, (J) increases with p\ up to a certain value and 
then gets saturated. This implies that below a certain 
strength, the defect strength is incapable of affecting the 
current. This picture changes dramatically when large 
a and small (3 are taken into account. In this case, (J) 
increases with p\ up to an /-dependent value and then 



starts decreasing. The maximum current sustained by 
the system is considerable. While in the case a = 0.8 (3 = 
0.05, the current for the normal ASEP is (3(1-/3) - 0.05; 
here we observe that disorder can remarkably enhance 
the current to almost a doubled value around 0.12. 

This type of constructive behaviour can be explained 
on the same grounds as for figure (8). Qualitatively, the 
impurities do not allow the overflow of particles into the 
system bulk which otherwise would have led to conges- 
tion and current reduction. If pi is small, the strength 
of defects is sufficient to block the the inflow of particles 
and reduces the current. If p\ is high enough, large inflow 
a will dominate and (J) is reduced. At an intermediate 
Pi we have a maximal current. Analogous to low a high 
f3, when both a and (3 are large, increasing p\ leads to 
current increments. If the density of defects is high, the 
current's increase would be linear in p\. For small /, the 
increase in current is rather linear for small p\ and after- 
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wards becomes more smooth. To deal in some depth, we 
next sketch the dependence of (p) on the defect strength 
Pi. For small a and large (3 (see figure (16)), one interest- 
ingly observes that if the defect concentration is relatively 
small i.e., less than 0.02, the effect of decreasing the de- 
fect strength (increasing p\) is to reduce the density as 
intuitively expected. Based on our previous arguments, 
defects are more influential if their concentration is rel- 
atively small [23| . Therefore, in the small concentration 
regime, the weakening of the defects leads to a sharp de- 
crease in the density. Beyond a certain pi, the further 
weakening of defects does not affect (p). For defect con- 
centration / above 0.02, the behaviour of (p) undergoes 
a qualitative change. As observed in figure (16), (p) in- 
creases up to a maximum and then starts diminishing. 
The maximum value of (p) depends on / and ranges be- 
tween 0.11 and 0.25. The reason is twofold. First, for 
small input rate a and intermediate concentration of de- 
fects, strong defects are still capable of forming rather 
large high density regions behind them which results in 
high (p). The second reason is due to the blocking of 
the particles outflow. Although (3 is high, strong enough 
defects are able to reduce this high outflow rate and ef- 
fectively reduce it. Consequently, the bulk density rises 
up throughout the bulk. Below a certain strength, the 
defective sites, although their numbers are not so small, 
are not only incapable of forming high density regions 
behind them but also incapable of effectively reducing 
the output rate. Therefore, (p) becomes decreasing. We 
now consider the case where a is large but (3 is small 
(see fig. 17). Here the overall effect of decreasing the 
defect strength is enhancement of (p). When the defect 
strength is reduced, the particles can more easily enter 
the chain and this leads to an increase in the bulk den- 
sity. In the limit of weak strength p\ — ► 1, we recover the 
normal value p — 1 — (3. When both a, and (3 are large 
corresponding to the MC phase in the normal ASEP, we 
still observe that the effect of a reduction of the defect 
strength is to enhance the density. Defects can drasti- 
cally reduce the density if their concentration and their 
strength are both large. Otherwise, their influence is a 
slight reduction of the density. We note the type of den- 
sity increment is rather similar in figures (17) and (18). 



IV. UNIFORM DISTRIBUTION FUNCTION 

So far our investigation has been restricted to the case 
where the disorder strength was limited to only two val- 
ues. In order to obtain a complimentary insight into the 
nature of disorder effect, it would be noteworthy to con- 
sider the case where the defect's strength can be chosen 
from a continuous interval. For this purpose, we con- 
sider the uniform distribution function for the strength 
of defect. Here one has two parameters namely a and 
6 which are the first and last points of the distribution 
interval. Like the binary distribution, one can introduce 
two classes. In the first class, (p) = ^4^- — 1 which corre- 
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FIG. 11: Figure 19 : Average current (J) versus a for low a 
high /3, high a low and high a high f3. L = 200. 
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FIG. 12: Figure 20 : Average density versus a for low a high 
/3, high a low f3 and high a high f3. L = 200. 



sponds to the case having fast hopping sites greater than 
unity. In the second class, 6=1 while a (a < 1) de- 
termines the lower limit of the defect's value. Note that 
in the second class (p) = =4^ < 1. We now exhibit the 
results for the latter case i.e., 6=1. The following dia- 
gram depicts the dependence of (J) on the lower end of 
the interval a. All diagrams which are shown next have 
been obtained by Monte Carlo simulations. The number 
of disordered samples over which the averaging have been 
performed is 1000, and the system size is 200. 

According to figure (19), in low a-high (3 and high al- 
low (3 regimes, beyond a = 0.1, the impurities do not 
affect the current and the system can maintain the nor- 
mal ASEP values a(l — a) and (3(1 — (3) respectively. In 
contrast, for both a and [3 larger than 0.5, (J) is a smooth 
increasing function of a. The current reaches its normal 
ASEP value 0.25 in the limit a — > 1. The behaviour of 
the average bulk density (p) on a is shown in figure (20). 

Contrary to the current diagrams, here the density ap- 
proaches the normal ASEP value which depends on a and 
(3. The high a-high f3 regime has the weakest dependence 
on a and beyond a = 0.1, (p) will be independent of a. 
Figure (21) depicts the phase diagram in the case where 
the average rate of hopping is unity: (p) = = 1 for 
various values of a. Analogous to the binary distribution, 
the overall effect of disorder is to enlarge the size of the 
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FIG. 13: Figure 21 : Phase diagram of the disordered ASEP 
for a uniform distribution of hopping rates. The distribution 
characteristics are specified in the figures. All distributions 
have unit mean but the variances are different. 
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FIG. 14: Figure 22 : Size of the low and the high density 
phases versus the variance for uniform and binary distribution 
functions (a = pi, f = 0.5 and (p) = 1). L = 200. 



MC phase and shrinkage of low and high density phases, 
respectively. 

Figure (22), exhibits the size dependence of the LD 
(HD) phase 5 in terms of the variance of the distribution 
functions for both uniform and binary distribution func- 
tions. The mean of the distribution functions is set to 
unity. 

For the uniform distribution, the size increment of the 
MC phase shows a more rapid dependence on the vari- 
ance of the distribution function in comparison to the 
binary distribution. The reason is due to the fact that in 



the uniform distribution, the frequency of small-hopping 
sites close to the lower limit of the distribution interval 
is more than those in the binary distribution. 



V. SUMMARY AND CONCLUDING REMARKS 

Let us now summarize what has been explored in this 
paper. We have investigated the statistical characteris- 
tics of the asymmetric simple exclusion process in the 
presence of spatially uncorrelated quenched disorder in 
the hopping rates via extensive simulation and numer- 
ics. Our findings cover two different distributions of hop- 
ping rates: binary and uniform. The conventional three- 
phase structure of the normal ASEP remains unchanged. 
Generically, the disorder affects the phase diagram by en- 
larging the maximal current phase, which in turn leads 
to squeezing the low and high density phases. This is 
accompanied by an overall decrease (increase) in the cur- 
rents (densities). We have managed to numerically solve 
the mean-field equations. Monte Carlo simulations are in 
support of the mean field solutions. In brief, the current 
exhibits a diminishing behaviour in terms of the defect's 
concentration in the chain when the input rate is small 
and the output rate is high. Analogously, it decreases 
when both the input and output rates are relatively high. 
Unexpectedly, in the case when the input rate is large and 
the output rate is small, the current shows an increas- 
ing dependence versus the defect's concentration. This 
demonstrates the nontrivial interplay of spatial sitewise 
disorder with the drive. We have also examined the prop- 
erties of the ASEP under uniformly distributed spatial 
disorder. Although the phase structure is similar to that 
of with a binary distribution, we have identified distinc- 
tive features between them. Our study has been limited 
to disorder distribution functions with finite second mo- 
ments. We expect to observe substantial different types 
of behaviours for those distributions having a long tail. 
Work along this line is in progress. 
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